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ABSTRACT 

The evolution of the multipolar structure of the magnetic field of isolated neutron 
stars is studied assuming the currents to be confined to the crust. We find that except 
for multipoles of very high order (l> 25) the evolution is similar to that of a dipole. 
Therefore no significant evolution is expected in pulse shape of isolated radio pulsars 
due to the evolution of the multipole structure of the magnetic field. 
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1 INTRODUCTION 

Strong multipole components of the magnetic field have long 
been thought to play an important role in the radio emis- 
sion from pulsars. Multipole fields have been invoked for the 
generation of electron positron pairs in the pulsar magneto- 
sphere. For example, the Ruderman & Sutherland (1975) 
model requires that the radius of curvature of the field lines 
near the stellar surface should be of the order of the stel- 
lar radius to sustain pair production in long period pulsars. 
This is much smaller than the expected radius of curvature 
1LJ ' of the dipole field. Such a small radius of curvature could be 
.^H | a signature of either an extremely offset dipole (Arons, 1998) 
, or of a multipolar field (Barnard & Arons, 1982). Further, 
» I ■ soft X-ray observations of pulsars show non-uniform surface 
temperatures which can be attributed to the presence of a 
quadmpolar field (Page & Sarmiento, 1996). 

Magnetic multipole structure at and near the polar cap is 
also thought to be responsible for the unique pulse profile of 
a pulsar (Vivekanand & Radhakrishnan 1980, Krolik 1991, 
Rankin & Rathnasree 1995). The recent estimates that there 
should be several tens of sparks populating the polar cap is 
also best explainable if multipole fields dictate the spark ge- 
ometry near the surface (Deshpande & Rankin 1998, Rankin 
& Deshpande 1998, Seiradakis 1998). Significant evolution 
in the structure of the magnetic field during the lifetime of 
a pulsar may therefore leave observable signatures. If the 
multipoles grow progressively weaker in comparison to the 
dipole then one can expect pulse profiles to simplify with 
age and vice versa. 

The evolution of the magnetic fields in neutron stars in 
general is still a relatively open question. During the last 
decade, two major alternative scenarios for the field evolu- 
tion have emerged. One of these assumes that the field of 
the neutron star permeates the whole star at birth, and its 
evolution is dictated by the interaction between superfluid 



vortices (carrying angular momentum) and superconduct- 
ing fluxoids (carrying magnetic flux) in the stellar interior. 
As the star spins down, the outgoing vortices may drag and 
expel the field from the interior leaving it to decay in the 
crust (Srinivasan et al. 1990). In a related model, plate tec- 
tonic motions driven by pulsar spindown drags the magnetic 
poles together, reducing the magnetic moment (Ruderman 
1991a,b,c). 

The other scenario assumes that most of the field is gener- 
ated in the outer crust after the birth of the neutron star 
(Blandford, Applegate & Hernquist 1983). The later evolu- 
tion of this field is governed entirely by the ohmic decay of 
currents in the crustal layers. The evolution of the dipole 
field carried by such currents has been investigated in some 
detail in the recent literature (Sang & Chanmugam 1987, 
Geppert & Urpin 1994, Urpin & Geppert 1995, 1996, Konar 
& Bhattacharya 1997, 1998). These studies include field evo- 
lution in isolated neutron stars as well as those accreting 
from their binary companions. The results show interesting 
agreements with observations lending some credence to the 
crustal picture. 

In this paper, we explore the ohmic evolution of higher order 
multipoles in isolated neutron stars assuming the currents 
to be originally confined in the crustal region. Our goal is 
to find whether there would be any observable effect on the 
pulse shape of radio emission from isolated pulsars as a result 
of this evolution. In section 2 we discuss the details of the 
computation and in section 3 we present our results and 
discuss the implications. 



2 COMPUTATIONS 

The evolution of the magnetic field, due to ohmic diffusion, 
is governed by the equation (Jackson 1975) : 
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™=-£ F Vx(IxVxB) I (1) 

where o(r,t) is the electrical conductivity of the medium. 
Following Wendell, Van Horn & Sargent (1987) we introduce 
a vector potential A = (0, 0, A,f,) assuming the field to be 
purely poloidal, such that: 

S{r,6,t) = -r sin6 Af(r,$,t), 

where S(r, 9, t) is the Stokes' stream function. S can be sep- 
arated in r and 9 in the form : 

S{r,9,t) = y^ y Ri(r,t) sin9 P{(cos9), 
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where Pj (cos(9)) is the associated Legendre polynomial of 
degree one and Ri is the multipole radial function. From 
equation (\w we obtain : 



d 2 Ri _ + _ AirRla dRi 

dx 2 x 2 1 c 2 dt 



(2) 



where x = r /R, is the fractional radius in terms of the stellar 
radius R* . The solution of this equation with the boundary 
conditions : 



OX X 

Ri = 



1 



at x — x c (3) 

for a particular value of I gives the time-evolution of the 
multipole of order I. Here, the first condition matches the 
correct multipole field in vacuum at the stellar surface and 
the second condition makes the field vanish at the core-crust 
boundary (where r — r c , the radius of the core) to keep the 
field confined to the crust. We assume that the field does 
not penetrate the core in the course of evolution, as the core 
is likely to be superconducting. 



2.1 Crustal Physics 

The rate of ohmic diffusion is determined mainly by the 
electrical conductivity of the crust. The conductivity of the 
solid crust is given by 

1 _ 1 1 

O <7 p h CTi mp 

where cr p h is the phonon scattering conductivity, which we 
obtain from Itoh el al. (1984) as a function of density and 
temperature, and the impurity scattering conductivity <Ji mp 
is obtained from the expressions given by Yakovlev & Urpin 
(1980). 

We construct the density profile of the neutron star in ques- 
tion using the equation of state of Wiringa, Fiks & Fab- 
rocini (1988) matched to Negele & Vautherin (1973) and 
Baym, Pethick & Sutherland (1971) for an assumed mass of 
1.4 Mq. As conductivity is a steeply increasing function of 
density and since the density in the crust spans eight orders 
of magnitude the conductivity changes sharply as a function 
of depth from the neutron star surface. Thus the deeper the 
location of the current distribution, the slower is the decay. 
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Figure 1. The evolution of the surface magnetic field for var- 
ious multipoles due to pure diffusion. The numbers next to the 
curves correspond to respective orders of multipole. All the curves 
correspond to Q = 0.0 and a depth of current concentration at 
Xo = 0.98 i.e., a density of p = 10 11 gem -3 . 



The thermal evolution of a neutron star has been com- 
puted by many authors, and it is clearly seen that the in- 
ner crust (p > 10 10 gcm~ 3 ) quickly attains an isothermal 
configuration after birth. At outer regions of the crust, the 
temperature follows an approximate relation, 



1/4 



T, p 



Pb 



(4) 



where T; is the temperature of the isothermal inner crust and 
Pb is the density above which the crust is practically isother- 
mal. As the star cools, larger fraction of the crust starts 
becoming isothermal, with pb being approximately given by, 



pb = 10 1 



VioV 



(5) 



The relations ^ and [| above have been obtained by fitting to 
the radial temperature profiles published by Gudmundsson, 
Pethick & Epstein (1983). For the time evolution of Tj we 
use the results of Urpin & van Riper (1993) for the case of 
standard cooling (the crustal temperature T m in their nota- 
tion corresponds to Ti above). 

A third parameter that should be considered in determin- 
ing conductivity is the impurity concentration. The effect of 
impurities on the conductivity is usually parametrised by a 
quantity Q, defined as Q — ^ — 1 where n is 

the total ion density, m is the density of impurity species 
i with charge Zi, and Z is the ionic charge in the pure 
lattice (Yakovlev & Urpin 1980). In the literature Q is as- 
sumed to lie in the range 0.0 - 0.1. But statistical analyses 
indicate that the magnetic field of isolated pulsars do not 
undergo significant decay during the radio pulsar life time 
(Bhattacharya et al. 1992, Hartman el al. 1997, Mukherjee 
& Kembhavi 1997) . It has been shown (Konar 1997) that to 
be consistent with this impurity values in excess of 0.01 are 
not allowed in the crustal model. 



Another important factor in determining the conductivity 
is the temperature of the crust. In absence of impurities 
the scattering of crustal electrons come entirely from the 
phonons in the lattice (Yakovlev & Urpin 1980) and the 
number density of phonons increases steeply with tempera- 
ture. The thermal evolution of the crust therefore plays an 
important role in the evolution of the magnetic field. 



2.2 Numerical Scheme 

To solve equation (Q) we assume the multipole radial pro- 
file used by Bhattacharya & Datta (1996, see also Konar 
& Bhattacharya 1997). This profile contains the depth and 
the width of the current configuration as input parameters 
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Figure 2. The ratio of the dipole surface field to the multipolc 
field is plotted as a function of age. The numbers next to the 
curves correspond to respective orders of multipole. All the curves 
correspond to Q = 0.0 and a depth of current concentration at 
x = 0.98 i.e., a density of p = 10 11 gem -3 . 



Figure 3. The ratio of the dipole surface field to that of the 
multipoles at 10 7 years as a function of Q. The numbers next to 
the curves correspond to respective orders of multipole. All curves 
correspond to a depth of x = 0.98 corresponding to a density of 
p = 10 11 gem -3 , at which the initial current is concentrated. 



and we vary them to check the sensitivity of the result to 
these. We solve equation (Q) numerically using the Crank- 
Nicholson method of differencing. We have modified the nu- 
merical code developed by Konar (1997) and used by Konar 
& Bhattacharya (1997) to compute the evolution of mul- 
tipolar magnetic fields satisfying the appropriate boundary 
conditions given by equation (0). 



3 RESULTS AND DISCUSSION 

In figures ^] and Q we plot the evolution of the various 
multipole components of the magnetic field, assuming the 
same initial strength for all, with time due to pure diffusion 
in an isolated neutron star. It is evident from the figures that 
except for very high multipole orders (Z> 25) the reduction 
in the field strength is very similar to that of the dipole com- 
ponent. For a multipole of order I there would be 2 reversals 
across the stellar surface. For typical spin-periods the size 
of the polar cap bounded by the base of the open field lines 
is ~ 0.01% of the total surface area. To contribute to the 
substructure of the pulse therefore the required multipoles 
must have a few reversals in the polar cap which demands 
that the multipole order must be five or more. On the other 
hand if the multipole order is very large (I > Z max ~ 20) 
the fine structure would be so small that it would be lost in 
the finite time resolution of observations. Therefore, I val- 
ues in the range 5 to Z max would be the major contributors 
to the observed structure of the pulse profile. However, as 
seen from figures [ jj and [0 multipoles of such orders evolve 
similarly to the dipole. Therefore no significant evolution 
is expected in the pulse shape due to the evolution of the 
multipole structure of the magnetic field. As discussed be- 
fore multipole orders contributing to the required field line 
curvature for pair-production are low, most prominently a 
quadrupole. As the evolution of these low orders are also 
very close to the dipole the radii of curvature of the field 
lines on the polar cap are not expected to change signifi- 
cantly in the lifetime of a radio pulsar. 

To test the sensitivity of these results on the impurity con- 
centration of the crust and the density at which the initial 
current is concentrated we have evolved models with vari- 
ous values of these parameters. The results are displayed in 
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Figure 4. The ratio of the dipole surface field to that of 
the multipoles at 10 7 years as a function of depth. The 
points marked in the plots here correspond to densities p = 
10 13 - 5 , 10 13 , 10 12 - 5 , 10 12 , 10 11 - 5 , 10 11 , 10 10 5 , 10 10 , 10 9 - 5 , 10 9 gem" 
The numbers next to the curves correspond to respective orders 
of multipole. All curves correspond to Q = 0. 



figures [g] and |4j where we plot the ratio of the dipole to 
higher multipoles at an age of 10 7 years. It is seen that the 
results are insensitive to these parameters, particularly for 
low orders of multipoles of interest. 

Krolik (1991) and Arons (1993) conjectured that except for 
multipoles of order l> i?*/Ar the decay rates would be sim- 
ilar due to the finite thickness Ar of the crust over which 
the current is confined. The evolution plotted in figure |lj 
assumes that Ar = 1.2 km for which R,/Ar ~ 8. However 
it is seen from figures |jj and Q| that significant decay oc- 
curs only for l>25, much greater than R t /Ar. This is most 
likely caused by steep increase in conductivity towards the 
interior. 



In conclusion, our results indicate that for a crustal model 
of the neutron star magnetic field there would be no signif- 
icant change in the multipolar structure with age. This fact 
seems to be corroborated by observations: studies identify- 
ing multiple components in pulse profiles (Kramer et al. , 
1994) show that the number of components does not vary 
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with the age of the pulsar. Thus the evolution of the multi- Wiringa R. B., Fiks V., Fabrocini A., 1988, Phys. Rev. C, 38, 

polar structure of the magnetic field is unlikely to leave any 1010 

observable signature on pulsar emission. This is in contrast Yakovlev D. G., Urpin V. A., 1980, SvA, 24, 303 
with the predictions from the plate-tectonics model of Ru- 
derman (1991a,b,c) which suggests a major change in the 
field structure with pulsar spin evolution. 
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